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On account of the Abel-Galois no-go theorem for the algebraic solution to quintic and higher 
order polynomials, the eigenvalue problem and the associated characteristic equation for a general 
noise dynamics in dimension d via the Choi-Jamiolkowski approach cannot be solved in general 
^ , via radicals. We provide a way around this impasse by decomposing the Choi matrix into simpler, 

i!^ . not necessarily positive, Hermitian operators that are diagonalizable via radicals, which yield a 

' set of 'positive' and 'negative' Kraus operators. The price to pay is that the sufficient number of 

Kraus operators is d* instead of d'^, sufficient in the Kraus representation. We consider various 
applications of the formalism: the Kraus repesentation of the 2-qubit amplitude damping channel, 
the noise resulting from a 2-qubit system interacting dissipatively with a vacuum bath; defining the 
maximally dephasing and purely dephasing components of the channel in the new representation, 
, and studying their entanglement breaking and broadcast properties. 



PACS numbers: 03.65.Yz, 03.67.-a 



I. INTRODUCTION 



^ I Any practical use of a quantum operation involves taking into account the effect of the ambient environment, and 
' the systematic study of such an effect is called open quantum systems. This now prevades a vast arena of studies, see 
2P , for e.g., the Refs. [H, [3| for two distinct flavors of the subject. The effect of the environment, interchangeably called 
' the reservoir or the bath, effects the system dynamics, in general, in two ways depending upon the commutability of 
• , the system and interaction Hamiltonian. If the two commute, then the process is a quantum nondemolition one, that 
is, there is dephasing without any energy exchange while if they do not commute, then there is dephasing along 
'J [ with dissipation These effects have been brought within the ambit of practical implementation by a number of 
I . very impressive experiments, for e.g., Q involving ion traps and @ using high-Q cavity quantum electrodynamics. 

As a result, the use of ideas from open quantum systems have become widespread in quantum information processing 
. A very useful tool in this regard is the Kraus representation which encodes the effect of the environment on the 
system of interest. Since a quantum operation that can be represented in the Kraus representation is guaranteed to 
5^ , be completely positive (CP), it is of importance to find the Kraus representation pertaining to different open system 
• models. The main aim of this paper is to develop Kraus representations for general N-qubit open system models, and 
use it specifically on a two-qubit model much discussed in the literature Q . The model considered is that of two qubits 
interacting with a bath, for e.g., an electromagnetic field in a squeezed thermal state, via the dipole interaction, which 
has been considered in detail for both pure dephasing [lo| as well as dissipative [llj system-reservoir interactions. 
The system-reservoir coupling constant is dependent upon the position of the qubit, leading to interesting dynamical 
consequences. Basically this allows a classification of the dynamics into two regimes: the independent (or localized) 
decoherence regime, where the inter-qubit distances are such that each qubit sees an individual bath or the collective 
decoherence regime, where the qubits are close enough to justify a collective interaction with the bath. 

When we attempt to obtain the Kraus operators, for both the independent as well as collective regimes for a version 
of the dissipative model [T]| . described below, we encounter a problem which has its origin in the famous Abel-Galois 



'Electronic address: |omkar@poornaprajna.org| 
tElectronic address: srik@poornaprajna.org 
t Electronic address: subhashish@iitj .ac.in 



2 



irreducibility theorem. We are able to circumvent this mathematical obstruction, in the case considered by resorting to 
the inherent symmetries in the model. This is accomplished by introducing the concept of extended Kraus operators. 
We can thus coin the word Abel-Galois integrable for such systems. It would be of interest, here, to note that the 
label Abel-Galois integrable could be ascribed to infinite dimensional systems also, such as the dissipative harmonic 
oscillator [l^l, which due to their inherent symmetries can be solved by quadratures. 

The plan of the paper is as follows. In the next section, we briefly discuss, the Abel-Galois irreducibility theorem, 
due to its relevance to our work. We then, in anticipation of its need for the Kraus operators of the general two- 
qubit dissipative model, introduce the concept of extended Kraus operators. This is followed by some illustrative 
applications to a typical single-qubit channel. The use of this formalism is, of course, not needed for the single qubit 
case, where it is possible to obtain the usual Kraus operators, but is intended to serve as an illustrative example to 
the main application of extended Kraus operators to two-qubit dissipative interaction with a vacuum bath, the 2AD 
(amplitude damping) channel. We then indicate why this would not work for the case of a bath at finite temperature 
and bath squeezing, the two-qubit squeezed generalized amplitude damping channel (2SGAD). Next, we discuss some 
features of the 2AD channel. Finally, we make our conclusions. 

II. ABEL-GALOIS IRREDUCIBILITY THEOREM 

The first famous no-go result in algebra, the irreducibility theorem, discovered independently by mathematicians 
Niels Henrik Abel and Jean Evariste Gallois (and anticipated earlier by Ruffini) states that polynomial equations 
of degree 5 or higher do not in general have solutions that can be expressed algebraically, i.e., in terms of addition, 
subtraction, multiplication, division and taking roots to a given finite order over the equation's coefficients and rational 
numbers (and any finite number of irrationals). Obviously there are infinitely many examples where solutions thus 
expressible do exist, a trivial example being /^{x) = {x — c)^ = 0, where c is a real or complex number, and where 
the solution x = c is manifest. Similarly solutions to a product of a quartic and a linear polynomial, f4{x)fi{x) = 0, 
can be obtained by solving those polynomials separately. For the general case of n > 5, one would have to resort to 
numerical methods like the Laguerre method or the Newton-Raphson method. 

Briefly, the argument, which is part of Galois theory, runs thus [l3| : a polynomial equation over rational numbers 
(or more generally, over the base field of given constants) admits a solution by radicals precisely if its Galois group 
is a solvable group. For polynomials upto quartic degree, the associated Galois group is solvable, but for quintic and 
above, unsolvable cases exist. 

Let the base field be the set Q of all rationals (augmented by at most a finite number of irrational constants, which 
we ignore for simplicity). Suppose we are given a polynomial 

5 

f{x)=J2»r^' (^QM- (1) 

If Xfc are the solutions to f{x) = 0, then: 

/(x) =nLi(.T-xfe) eE, (2) 

where the solutions Xk exist in general in the splitting field E, the extension field of Q, which is the smallest subfield 
of C containing the roots of f{x). 

Multiplying out Eq. ^ and comparing it with Eq. ([T} shows the a^'s to be elementary symmetric functions of the 
roots Xj-. 

Pa = ^XjiPs = ^XjXk;P2 ^ ^ XjXkXi, 

Pi = 2^ ■ ■,Po = UjXj, (3) 

j ■' 

where aj = (— l)-'+^/3j. We observe that every permutation a of the roots, which is clearly an automorphism of E, 
leaves Pj e Q, and hence Q, fixed. It can be shown that this is the only automorphism of E that leaves Q fixed. 
The group of automorphisms ct of E such that a{q) = q, with g e Q is called the Galois group associated with the 
polynomial, and denoted Gal(E/Q) C S'„, the group of permutations on the roots of f{x). 

A subgroup N of group G is called normal if and only if gN = Ng for all g G G, oi equivalently, g^^Ng = N, 
that is the normal subgroup is invariant under conjugation of elements of G. The relationship is denoted N <iG. The 
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normal subgroup defines the quotient or factor group, denoted G/N. A subnormal series (or tower) of a group G is a 
sequence of subgroups: 

{e}<Ao<Ai<A2---< An = G, (4) 

terminating in the trivial subgroup. A composition series is a subnormal tower such that Aj ^ Aj+i and Aj is the 
maximal normal subgroup of Aj^i. The group G is solvable only if the factor groups Aj^i/Aj in the composition 
series are Abelian. The Abel-Galois theorem can be stated as follows: a polynomial equation is solvable by radicals 
iff its Galois group is solvable. 

For all n, the maximal normal subgroup of 5„ is An, the alternating group, which is the subgroup of even permu- 
tations on S. For n < 4, all subgroups of 5„ are solvable. However for n > 5, An is non-Abelian, as is the factor 
group An/{e}, so that Sn is not solvable, and we find that polynomials that are quintic or of higher degrees are in 
general not solvable. 

The significance of the Abel-Galois theorem for us is that our work on deriving the Kraus operators, for general 
two-qubit systems, requires diagonalizing a (density) matrix in a Hilbert space of dimension cP = 16, which involves 
solving the characteristic equation of a self-adjoint matrix of degree 16. In general, the relevant Galois group is 
Sie, which is unsolvable. We circumvent this problem, in any dimension, by making use of symmetry properties of 
the matrix to circumvent the problem of diagonalization. At times, it may happen that even if the Galois group is 
solvable, diagonalization can be so tedious that our alternative is preferable. 



III. EXTENDED KRAUS OPERATORS 



A transformation of a quantum state is a CP map if and only if it evolves a density operator according to the 
prescription p^ : 

(5) 

j 

where Aj are at most cP operators that satisfy the completeness condition = ^- Here we will study the 

conditions under which a CP map allows a description of the form: 

p ^ E{p) = p' = E A+p (A+)J - ^ A-p {A-)l , (6) 
j=i k=i 

where p + v > cP and the extended Kraus operators A^ must satisfy the new completeness condition 

j=i k=i 

Let 1^^) represent a vector 'unfolding' of A^ . Then the necessary and sufficient condition for the extended operator 
sum representation to be CP is that: 

E E IAT)(-4fe I ^B+-B-=B, (8) 

i=i fc=i 

where B is the Choi matrix for the map. 

To see this, we note that the Choi matrix corresponding to noise £ can be written as: 

B^Y.\^){k\®E{\j){k\). (9) 

3 k 

Let B — "Y^j Bj be any Hermitian partition oi B. The partition elements Bj can be spectrally decomposed with positive 
or negative eigenvalues. Let their corresponding eigenvectors, normalized to the absolute value of the eigenvalues be, 
1^^) and \A~), respectively. It follows from the properties of matrices that [l5| : 

\At){Af\ = J2 ® {A^\j){k\Af^)) . (10) 

3,k 
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Inserting these values in Eq. ([8]) and comparing the result with the expression in Eq. ([9]), we have: 



fe=i 



= J2\j){k\^£{\j){k\)., 

jk 



J2At (\j){k\AP 



(\j){k\Ar^ 



(11) 



from which Eq. © follows, because the sum- difference operation on each element |j)(fc| reproduces the effect of £ on 
that element, and thus also on any linear combination, such as p, of those elements. 

We consider an example of the above idea, applied to the generalized amplitude damp- 
ing channel (GAD) The Choi matrix corresponding to GAD channel with elements 



IS 



B = 



fl-X+pX \/l~A' 

pX 

(l-p)A 




\ Vl~A 
A decomposition of the above matrix is B = B^ — B^ , with 



(12) 



1-pXj 



B^ 



fl-X+pX + 



\ 4 





pX 

(l-p)A 

1 - pA -f 2 



l-A 



■,B- 



1^ 


2 






































V- 


4 








2 / 



(13) 



Both the matrices are seen to be Hermitian. The positive Kraus operators are 



\/4 + 2\/l- A 


- 2A- 




2 


•\/4 + 2Vl- A 


- 2A + 





2A(l-2p) + v^ 
3\/T^ 





y/{l-p)X 





k: = 





0\ 
pA j ' 



(14) 



where a = 9(1 — A) + 4A^(1 — 2p)^. The negative Kraus operators are 



i ri ' 



^2^ = 



V3{1 - A)i/4 



(15) 



Of course, the extended Kraus formalism used here is not necessary, since the usual Kraus operators for these 
channels can be obtained analytically. The purpose here was to serve as a simple example that sets the scene, in what 
follows, for the non-trivial application of the formalism in the two-qubit case. 



IV. INTRODUCTION TO TWO-QUBIT DYNAMICS 

Consider the Hamiltonian, describing the dissipative interaction of two qubits with the bath via the dipole interaction 



as 



i 



H — Hs + Hr + HsR 

2 



E f^nS^r, + E ^k{h\bus + 1/2) - ^^Y.Y.^^^^■hsi^^)ist + s-)h^^ - h.c:_ 



(16) 
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Here Hs is the system, Hn the bath, and Hsr the interaction Hamiltonians, respectively, and /2„ are the transition 
dipole moments, dependent on the different atomic positions r„. Also, are the dipole raising and lowering 
operators, respectively while is the energy operator of the nth atom, and bl^, b-^^ are the creation and annihilation 

operators of the (bath) field mode ks with the wave vector fc, frequency ujk and polarization index s = 1, 2 with the 
system-reservoir (S-R) coupling constant 

9ks 



hArn) = ( 



2ehV 



ik.rn 



(17) 



In Eq. (jl7p V is the normalization volume and e^^ is the unit polarization vector of the field. It can be seen from Eq. 
(|17p that the S-R coupling constant is dependent on the atomic position r„ leading to the possibility of considering 
the dynamics in the independent or collective regimes, depending on whether the qubits are far apart or close with 
respect to the length scales set by the environment. Assuming separable initial conditions, and taking a trace over 
the bath the reduced density matrix of the two-qubit system can be obtained [ll| . We will now attempt to obtain 
the Kraus operators for this model, in a unified way for both the independent as well as collective regimes, first for 
the case of a vacuum bath, that is, the 2 AD channel and then, point out the difficulties when we encounter finite 
temperature and bath squeezing. 

The two-qubit density matrix in dressed state basis is 



P = 



( Pee Pes Pea Peg 

Pse Pss Psa Psg 

Pae Pas Paa Pag 

\ Pge Pgs Pga Pgg 



(18) 



The time-evolution of the two-qubit density matrix p to p' = £{p) is given by: 



Apee J Pes M p^a 

J* Pse Bpss + Cpee PPsa 

M*pae P*Pas D Paa + Epee 

L*pg, T*pgs + {U* -iV*)pse Q*Pga + {~lS* 

where A, B,C, . . . ,U, V, S, R are given Appendix [X] 

The Choi matrix for the above interaction, which is the density operator {£ (g) I)(|<I>) (<I>|), where |$) = |00)|00) 
|01)|01) + |10)|10) + |11)|11), is given by: 



Lpeg 

Tpsg + {U + lV)pe 
QPag + {iS - R)p e 
)Pae Pgg + Fpss + Gpaa 



Hpee / 



(19) 



B = 



( ^ 
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U + iV 
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-iS* - R* 
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\ L* 




























Q* 














1 / 



(20) 



It turns out (as can be found using Mathematica software) that the relevant characteristic equation is cubic, so that 
the Galois group is a subset of S3. This is so because the matrix in Eq. (1^ is sparse, and is found not to be the 
case otherwise. However the solutions are so tediously long that they would be hardly of practical interest, and it is 
advantageous to use our method, which is applicable quite generally (even for non-sparse Hermitian matrices). 

As a result, the formalism of Section (jllip will be used to derive extended Kraus operators. We find the following 
decomposition convenient: 



B = Z?diag + Bj + Bm + Bl + Bp + Bq + Bt + Bu + Bv + Br + Bs, 



(21) 



where Sdiag is the submatrix of B consisting of precisely the diagonal entries in B, and O's elsewhere; Bj is the 
submatrix consisting only of the conjugate terms J, J*, and O's elsewhere, and so on. 
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A. Diagonal terms 

It is straightforward to see that by diagonahzing Bdiag = (A|0000)(0000| + C|0001) (0001| + £:|0010)(0010| 
if |0011)(0011| + S|0101)(0101| + F|0111)(0111| + DjlOlO) (10101 + G|1011)(1011| + |1111)(1111|), we have 



diag 



(22) 



Then the Kraus operators obtained by 'folding' each eigenvector /C^ is: 



/O ON 





Vl Oy 

/O ON 

10 

Vo 0/ 

a 0' 



vO 



K 



D 



^0 ON 





VO 1 0/ 

/O ON 



10 

Vo 0/ 



c 



^0 o^ 





vo 1 oy 

^0 ON 

10 



vO Oy 



(23) 

















<l 


















;) 




\l 











^0 ON 

10 



vO Oy 



These operators lead to the evolution £diag, which transforms only the diagonal elements, killing off the rest: 



Apee 

+ ^^t _ I Bpss+Cpee 







Dpaa + Epee 



(24) 



Pgg + Epss + Gpaa + Hpe 



j € T = {H, G, F, E, D, C, A, B, 1}. 



B. Off-diagonal terms 

By diagonalizing Bj = J|0000) (0100| + J*|0100)(0000|, we have: 

6,7^|/c+)(/c+|-|/C7)(/C7| 

By 'folding' the vectors, we obtain the Kraus operators 













\J\ 


-e"*J 








2 













Vo 









(\ ON 

e'"^' 



\0 0/ 



(25) 



(26) 



where = —(ujq — ^V2)t- The evolution produced by these operators transforms two conjugate elements of p into 
the two corresponding elements of p', while annihilating all other elements in the density operator. Thus: 



Pj 



KjP{K-,)^ 



/ e 

e-'^-'\J\p 


V 



"f-'\J\pes 0\ 





0/ 



(27) 
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Proceeding thus with the other terms in Eq. (|2ip . we obtain pairs of positive and negative Kraus operators: 



/I ON 

|M| [ 

2 ±e">"'' 

\0 0/ 



/O 
1 








VO ie*-* 




/O 

1 



\0 ie'-^^y 




10 





±6^-^ 







\0 ±e"^^ 




/O ON 
10 

±e"^p 

Vo oy 

0\ 



o/ 




0/ 



\0 ±: 



(28) 



where = -2woi, 0m = -(wo + ^^12)^, 0p = -2rii2i, 0t 

-(wo - ^12)t, 4>S = -(wo - ^12)t + 7r/2, 0ij = -(wo - fJi2)t + TT. 

These operators lead, analogously to Eq. (|27)) to the partial evolutions p^j, p^, pp, p^, p^, / 



= Pdiag + + + p'l + p'p + p't + P'u + P'v + p'q + p'r + p's^ 



{LU0 + ni2)t + TT/2,(l>Q = 

<'q and py, which satisfy: 
(29) 



while 



B = ^|/C+)(/C+| + ^(|/C+)(/C+|-|/C-)(/C7|), 



(30) 



where S = {J, M, X, P, T, [/, Q, i?, S}. 

In short, our strategy to circumvent the Abel-Galois theorem is to replace the problem of diagonalizing the Choi 
matrix, by that of diagonalizing simpler Hcrmitian matrices that sum up to the Choi matrix. The elements of this 
Hermitian partition are diagonal or rank-two matrices, and thus trivially or readily diagonalized, in the case we 
considered. The price to pay is that the number of Kraus operators required to implement the CP map on a density 



operator in an d-dimensional Hilbert space can be as many as — (P 



general more than d? , the sufficient number of Kraus operators in the conventional formalism Q . 



2 (cf. Section irVC|) . which is in 



C. Abel-Galois non-integrability 

In our method, Sdiag is a diagonal matrix while Bi {i = J, M, L, P, Q, T, U, V, i?, S) are all rank-2 Hermitian matrices, 
which are readily diagonalizable. We will assume that the base field is Q{B), i.e., Q augmented by the entries in B. 
In that case, the Galois group of Bdmg is trivial (consisting of the identity permutation) in that the splitting field of 
the characteristic equation 

Xdiag = ix-A)ix- C)ix - E){x -H){x- B){x - F){x - D){x ~ G){x ~ 1), (31) 

is the base field itself. The ten remaining component matrices of B are rank-2 matrices. The characteristic equation, 
for example, for Bj in Eq. (psjl . is: 

X., = x^~\J\\ (32) 

with solutions ±| J| = zLy^ j'^ + jj , which are in general irrational. Hence the symmetry group is 5*2 consisting of the 
identity element and an interchange. 

Thus our method of circumventing the Abel-Galois theorem can be considered as reducing a problem requiring the 
solution of (the unsolvable) Sj_2 to one requiring that of {82)^'^ ^'^ j|- ^^j-jjg thai for the particular form pp]) 

of the Choi matrix for the 2AD channel, the characteristic equation is cubic and indeed solvable (as can be found using 
Mathematica) . However, the resulting eigenvalues and eigenvectors are so ponderous, that our method is preferable. 

For the two-qubit squeezed generahzed amplitude damping (2SGAD) channel [Tlj, the problem is seen to be Abel- 
Galois non-integrahle, i.e., it does not admit analytic Kraus operators. 
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V. AN APPLICATION OF THE OPERATOR SUM-DIFFERENCE FORMALISM 

We indicate some features of the operator sum-difference formalism in identifying special cases of a given noise that 
may have special interest. 

A. The maximally dephasing component 



It may be interesting to note that the Kraus operators , j e T = {H, G, F, E, D, C, A, B, 1}, which correspond to 
the diagonal terms in the Choi matrix, by themselves constitute a CP trace-preserving (CPT) dynamics. Physically, 
this channel represents a dynamics of the populations (diagonal terms of the density operator) that is completely 
dephasing. The evolution under this channel, which may be called the maximally dephasing component (MDC) 
channel, is given by Eq. (|24|) . Because all Kraus operators of this channel are of rank 1, it has the property of being 
entanglement breaking. In any dimension, it has an analytical operator sum (as against, sum-difference) representation. 

A quantum channel £ is called entanglement breaking and trace-preserving (EBT) if given any input state F, the 
state (I(8)f)r is separable. This is equivalent to the separability of the Choi matrix {I ^ £)\ip'^){tp'^\ and also to the 
condition that it can be expressed in the Holevo form 



fBs(p) =^i?^Tr(^^,p), 



(33) 



where {Ri} is a set of density operators and {Fi} is a set of positive operator valued measures (POVM) (Tq . 

The 2 AD channel is asymptotically {t — > oo) EBT because only the Kraus operators elements of the MDC 
survive. In the asymptotic limit the channel is characterized by the fact that in I^U^, H,F,G — > 1, whereas 
A, B, C, D, E, J, K, L, M, P, Q, X, Y — !• 0. Therefore, only the following 4 rank-1 Kraus operators, which are elements 
of the MDC channel, survive: 



X+(oo) 



/OOOO 



\ 1 



i^+(oo) 



/o o^ 





Vo 1 oy 



/o o^ 





Vo 1 oy 



i^+(oo) 



/o o^ 





Vo 1/ 



(34) 



This channel produces the asymptotic state |5'o 



1 11), so that the corresponding Choi matrix has the separable 



form |11)(11|(8)/ and the channel is entanglement breaking: for an arbitrary initial state of the 2-qubit system possibly 
entangled with any other system, asymptotically the 2-qubit system factors out to |11). 

Furthermore, asymptotically it is of an extreme point of entanglement breaking, whereby the channel maps all 
input states to a single point in state space, the pure state |11): it is thus a point channel |16| . The Kraus operators 
of EB channel can be expressed as Ki = \/Ri\jk){lm\^/Fi, which are seen to reproduce the expression in Eq. (|33p . 
For the considered point channel R, ^ R ^ F, = {|00)(00|, |01)(01|, |10)(10|, |11)(11|} and the resulting 

are easily seen to conicide with the operators in Eq. 

A special case of the entanglement breaking channel is the quantum-to-classical (QC) measurement map, wherein 
the states Ri in Eq. (p3|) are orthogonal projectors \ej){ej \ of a fixed basis. Acting on any 2-system state, it produces 
a QC state, which has the form: 



(35) 



A fundamental result here, which refines the equivalent result for EBT channels, is that a channel A is of QC-type 
(meaning that (I (E" A)(pab) is a QC state for any bipartite state pab) if and only if the correponding Choi matrix 
is a QC state [l4|- A particular feature of interest here is that, for any QC-type channel A, given any orthonormal 
basis there exists at least one state p*i(t>) diagonal in this basis, which is A'^-copy spectrum-broadcastable using 

the channel. If is identified with the channel basis {|ej)}, then p*{e) is fully broadcastable. Here the broadcast is 
state fj shared between Alice and A^ other parties, such that the reduced density matrix at each party has the same 
eigenvalue spectrum as p*{(l>) (spectrum broadcast) or is identical with p*{(j)) (full broadcast). 

The question then arises of whether the MDC channel is also of QC-type. The answer, according the above quoted 
result, is in the negative, as seen from the matrix Sdiagj formed by the diagonal terms of B in Eq. (PO]) . which cannot 
be cast in the form (|33p with Ri given by orthogonal projectors. To see this, we note that because only diagonal 
terms are present, both aj and \ej) in Eq. (1351) must be diagonal in the computational basis. Then it is clear why 
^diag, though separable, is not QC. In particular, looking at the first 4 diagonal terms, we should have C = H = 0, 
which is not in general guaranteed. 
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B. The purely dephasing component 

By contrast, a non-dissipative and purely dephasing operation is obtained from a channel comprising the extended 
Kraus operators j G S = {J, M, L, P, T, U, V, Q, R, S}. This has the effect of dephasing the off-diagonal terms, but 
killing off the diagonal terms. Restoration of the diagonal components is obtained by use of positive Kraus operators 
given by all the projetors: Hoo = |00)(00|, Hoi = |01)(01|, Hio = |10)(10| and Hn = which collectively form 

the set U. The elements of this set are just the Kraus operators K^, K\, and in Eq. (|24|) . with D,A,B = 1. 

Together the operators in the set SUU constitute a CP map that is a purely dephasing component (PDC) of the 2AD 
channel, which leaves the populations (diagonal components) unchanged, but otherwise reproduces the (dephasing) 
effects of the 2AD channel. Its action is given by: 

ieu ies 

(Pee J Pes M Pea Lpeg 

J*Pse Pss PPsa Tp^g + {U + iV)pes 

M*pae P*Pas Paa QPag + (iS - R)pea 

L*Pge T*pgs + {U* ~iV*)p,e Q*Pga + {~iS* ~ R*)pae Pgg 

Unlike MDC, this is not entanglement breaking at finite time. To see this, suppose that the input state is V'in = 
^(|00, 00) + |11, 11)), i.e., a 1-bit entanglement between the given 2-qubit system, with another 2-qubit system. The 
effect of the MDC channel is seen to be 

Pout = ^ (|00, 00)(00, 00| + 111, 11)(11, 111 + Jie'^"- 100, 00)(11, 11| -I- e"*"^^ |11, 11)(00, 00|)) . 

The above state lives in four dimensional Hilbert space, given by H2 ^ "^2, spanned by kets {|00), |11)} in each 
two-qubit system. The state is thus entirely equivalent to a correlated state of two qubits given by: 

Pout = I (|0,0)(0,0| + |1,1)(1,1| + J(e'^-|0,0)(l,l| +e-^^-|l,l)(0,0|)) , 

for which the concurrence [l8l | , a measure of entanglement for a mixed 2-qubit system is J, implying that entanglement 
is not broken except asymptotically (when J — )■ 0) . 

VI. CONCLUSIONS 

In the problem of deriving the Kraus representation of a given noise dynamics via the Choi- Jamiolkowski approach, 
we have developed a method that circumvents the impasse due to the Abel-Galois no-go theorem for the algebraic 
solution to quintic and higher order polynomials. Our idea is to obtain a Hermitian decomposition of the Choi matrix, 
which yield a set of 'positive' and 'negative' Kraus operators, satisfying the completeness condition. 

The price to pay is that, in general, the sufficient number of Kraus operators for a d-dimensional system is c?^, 
rather than the usually sufficient number of d?. 

We have applied this formalism to three problems: the sum-difference representation of the 2-qubit amplitude 
damping channel; defining two mathematically interesting limits of any channel in this representation, namely the 
MDC and PDC channels, whose entanglement breaking property was studied. In particular, the MDC channel is 
EBT, whereas the PDC channel is not, except asymptotically. 

Appendix A: The two-qubit amplitude damping channel 

Information useful for the description of the 2AD channel, discussed in Sec. IV, are: 

A =. e-2rt: b = e-(r+ri2)t. c = ^±^(1 - e-(r-ri2)t)e-(r+ri2)t. jj = e-(r-ri2)t. e = fe^(l - 

g-(r+ri2)i)g-(r-ri2)t. p ^ i _ g-(r+ri2)i. q ^ i _ g-(r-ri2)f. j ^ g^z(ajo-Oi2)tg{3r+ri2)i/2. ^ ^-i2u;ot^~Tt.^ 
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M = e-^("o+^i2)tg-(3r-ri2)t/2. p ^ g-i2ai2tg-rt. q ^ g-i(wo-ai2)tg-(r-ri2)«/2. ^ g-i(wo+f2i2)tg-(r+ri2)t/2. 



+ 

R = 

S ^ 
U = 
V = 



£+£12 
2r 

r + ri2 
r-ri2 



r-ri2 



£+lH(i_g-(™) + £^ 



-(r+ri2)t 



r2 


+ 4r!?2 


r 


-ri2 


r2 


+ Anl^ 


r 


+ ri2 


r2 


+ 4r!?2 


r 


+ ri2 



r2 + Anj^ 

The terms wqi ^12, T, ri2 are as defined in [ill . 



(1 - e-(r-ri2)t) _ ^iii(l - g-^r*) 

g-.Kc.o-ni2)tg-(r^ri2)t/2 [2f]i2e~" sin(2f]i2t) + r(l - e"" cos(217i2i))] , 
g-i(c.o-ni2)tg-(r-ri2)t/2 [2ni2{l - e"" cos(217i2t)) - Te"" sin(21]i2t))] , 
g-^(^o+Ol2)^g-(^+^l2)t/2 [21^126-" sin(2r!i2t) + r(l - e"" cos(2r!i2i))] , 
g-i(c.o+ni2)tg-(r+ri2)t/2 [21^12(1 - e"" cos(217i2i)) - Te-" sin(21]i2t))] . 
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